On the relevance of center vortices to QCD 
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In a numerical experiment, we remove center vortices from 
an ensemble of lattice SU(2) gauge configurations. This re- 
moval adds short-range disorder. Nevertheless, we observe 
long-range order in the modified ensemble: confinement is 
lost and chiral symmetry is restored (together with trivial 
topology) , proving that center vortices are responsible for both 
phenomena. As for the Abelian monopoles, they survive but 
their percolation properties are lost. 

PACS numbers: ll.15.Ha, 12.38.Aw, 12.38.Gc, 11.30.Rd 



The essential non-perturbative properties of QCD are 
confinement and chiral symmetry breaking (xSB). It has 
been observed through numerical lattice simulations that 
these two properties persist in the quenched theory up to 
a critical temperature T c ~ 22QMeV, where confinement 
is lost; chiral symmetry appears to be simultaneously re- 
stored 0. The disorder which leads to the area law for 
the Wilson loop thus seems to be somehow tied to the ex- 
istence of a chiral condensate. Although effective mecha- 
nisms have been proposed to explain confinement or xSB, 
no successful common explanation is yet available. 

Two effective descriptions of QCD have been receiv- 
ing a lot of attention: one considers instantons as the 
effective degrees of freedom (d.o.f.), the other chromo- 
magnetic monopoles. Instantons are natural candidates 
to explain xSB: each instanton is associated with a zero 
mode of the Dirac operator [|j , and there must be an ac- 
cumulation of zero eigenvalues to obtain a quark conden- 
sate Q. Above T c , the vanishing of the condensate must 
correspond to a qualitative change in the instanton en- 
semble (see, e.g., ||). On the other hand, it is unlikely 
that instantons play a significant role in confinement (see 
|^| for a recent discussion). An attractive mechanism for 
confinement is dual superconductivity of the QCD vac- 
uum ||. Considerable evidence for this dual Meissner 
effect has been accumulated on the lattice, including a 
disorder parameter demonstrating the condensation of 
chromo-magnetic monopoles below T c This conden- 
sation has been observed directly after gauge-fixing to 
"Maximal Abelian Gauge" ||, and the idea of "Abelian 
dominance" has emerged, according to which the Abelian 
d.o.f. of the Yang-Mills field encode all its long-distance 
(IR) properties. Indeed, %SB and its restoration have 
been observed in the Abelian sector 0. The Abelian 
dominance scenario has some flaws, however: it does 
not explain the breaking of the adjoint string, and the 



Abelian string tension differs slightly from the Yang-Mills 
one Moreover, d.o.f. more elementary than Abelian 
monopoles, embedded in them and solely responsible for 
the physics assigned to them, cannot be ruled out. 

The idea of center vortices, which initially failed due to 
their misidentification jljj ], has been successfully revived, 
first as an embedded model inside the Abelian sector [jl2| , 
then without reference to Abelian projection [13). Cen- 
ter vortices are exposed by gauge-fixing: after a gauge 
transformation which brings each lattice link as close as 
possible to a center element of the gauge group, vortices 
consist of defects in the center-projected gauge field. The 
idea of center dominance is again that the center d.o.f. 
encode all the IR physics. The density of center vortices 
seems to be a well-defined continuum quantity |h^ . |l4f : 
the center string tension appears to more or less match 
the original one; and an explanation for the behaviour of 
the adjoint potential has been proposed |jl5). However, 
chiral symmetry has not yet been studied in this context. 

An additional problem with Abelian and center domi- 
nance is that the relevant d.o.f. are identified only after 
gauge-fixing. Gauge-fixing non-Abelian fields is notori- 
ously ambiguous, and different Gribov copies produce dif- 
ferent Abelian monopoles or center vortices, whose prop- 
erties like the string tension differ slightly. For this rea- 
son, we are a priori suspicious of effective models which 
involve gauge-fixing, and so we designed a simple nu- 
merical experiment to disprove the center-dominance sce- 
nario. For simplicity, we consider the gauge group SU(2), 
with center i?2- Our starting point is an ensemble of lat- 
tice gauge fields representative of the continuum. We 
identify center vortices in this ensemble, and construct 
from it a modified ensemble where all center vortices have 
been removed by flipping the sign of a subset of SU{2) 
gauge links. This operation introduces a lot of disorder 
in the gauge field. Nonetheless, these disordered gauge 
fields now have a trivial, vortex-free center projection, 
and so according to the credo of center dominance they 
should not confine. That is, our introduction of short- 
range disorder should at the same time bring long-range 
order. To our surprise, this is indeed what happens. One 
may then ask if the spectral properties of the Dirac op- 
erator are not also dominated by the center components 
of the gauge field. In that case, our modified ensemble 
should show no sign of xSB, since its center projection is 
the trivial (perturbative) vacuum. Indeed, this is what 
we observe: removal of center vortices causes both loss of 
confinement and restoration of chiral symmetry. 
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The next intriguing question regards the fate of 
Abelian monopoles as center vortices are removed. They 
do not disappear; on the contrary, the introduction of 
short-range disorder increases their number. However, 
we observe the complete disappearance of monopole cur- 
rent loops winding around the periodic lattice: we can 
thus identify these as the fundamental objects associated 
with confinement in the Abelian sector, apparently influ- 
enced by the underlying center d.o.f. 

Finally, we investigated the effect of Gribov copies 
which caused our initial skepticism. We repeated our 
experiment on the same SU(2) ensemble, but introduced 
a systematic tolerance in the gauge condition to be satis- 
fied before identifying the center vortices to be removed. 
Although the location and number of center vortices we 
removed varied appreciably, the modified SU (2) ensem- 
ble was always non-confining and chirally symmetric. 

The numerical experiment - We start from an ensemble 
of SU(2) lattice gauge fields obtained by Monte Carlo 
using the standard Wilson plaquette action. To identify 
center vortices, we gauge-fix our configurations in order 
to bring each SU(2) gauge link U^(x) as close as possible 
to an element of the center Z-i = {+1,-1}. We therefore 
try to iteratively maximize 



Q({u,})^j2( TTU ^y 



(i) 



as in [|13|, where this gauge is called the "direct maxi- 
mal center gauge." The gauge-fixed SU(2) links, denoted 
Uf* (x), are then projected to Z2 elements Z^{x) using 



Z„{x) = sign[Tr U° F {x)\ . 



(2) 



Plaqucttcs in the ^-projected theory with value —1 rep- 
resent defects of the Zi gauge field called P-vorticcs [jl2] . 
Numerical evidence has been presented flj,|l3| show- 
ing that plaquette-like P-vortices signal the presence of 
macroscopic, physical excitations, called center vortices, 
in the unprojected original SU (2) configuration. 

Consider then the modified SU{2) configuration made 
of gauge links U'(x) constructed as 



K(x) 



Z^x) U,(x) 



(3) 



The gauge transformation which maximizes Q({U^}) 
in (P also gives the same maximum to Q({U'}), so 
that the modified gauge-fixed links U'^ F {x) are simply 
Zfj,(x) U^ F (x). Therefore, we instantly know the center 
projection Z' (x) of U'^{x): 

Z'^x) = sign[Tr U'f F (x)} = sign[Tr Z M (x) U^ F (x)} 

= (sign[Tr^ GF (x)]) 2 = +l. (4) 

Every modified configuration U' thus projects onto the 
trivial Zi vacuum: all center vortices have been removed. 



Our ensemble consists of about 1000 SU (2) configura- 
tions, on a 16 4 lattice at f3 = 2.4. To maximize Q({?7^}) 
(Eq.|l]) we use standard overrelaxation, stopping when 



e = J]A(Tr U^{x)f < 10" 



(5) 



from one gauge-fixing sweep to the next. 




FIG. 1. Normalized plaquette distribution on the original 
and modified ensembles (SU(2), f3 = 2.40). Center- vortex 
removal increases short-range disorder. 

In Fig.l we show the distribution of SU(2) plaquette 
values on the original and modified ensembles. It is ap- 
parent that under the sign flip Eq.(||), many SU(2) pla- 
quettes acquire a negative value. The modified ensemble 
has an increased action, i.e. more short-range disorder. 
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FIG. 2. Creutz ratios on the original and modified ensem- 
bles. The dashed band is the string tension result of [jhi|. 

Results - In Fig. 2 we present our results for the Creutz 
ratios xb„r = ~H( w R,r)( w R-i,R-i)/( w R,R-i) 2 } con- 
structed from averages (Wr,t) of R by T Wilson loops on 
the original and modified ensembles. For large R, XR,R 
tends to the string tension a. On the modified ensemble, 
the Creutz ratios clearly decrease and tend to zero. De- 
spite the increased short-range disorder, long-range order 
has been created and confinement has been lost. 



2 



10 



: ' ' ' ' I ' ' ' ' I ' ' ' ' f ' V ' t ' ' 
^Original ensemble + + + + + 



P5 6 '— K=a I 

f ,< '■ + 

t 4 — .•>■■:'.*: 



10 

A B 

e-j : 

CB 6 — 
^ : 



Modified ensemble 







2 



i 



6 



8 



10 



T 



FIG. 3. Wilson loop values (Wr,t) on the original and 
modified ensembles. Note the parallel lines for successive R 
in the latter: upon center-vortex removal, confinement is lost. 

This is even clearer if one looks directly at the Wilson 
loop values. In Fig. 3 we show — \&(Wr } t) as a function 
of T. For a fixed R, points at successively larger T form 
a line whose asymptotic slope is the value of the static 
potential V(R). The lines corresponding to the modified 
ensemble are parallel, indicating that V(R) does not grow 
with R: the string tension has vanished. 
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FIG. 4. Quark condensate (ifnp) (m q ) on the original and 
modified ensembles. The dashed line corresponds to "poor" 
center- vortex identification (see text). In all cases, cen- 
ter-vortex removal restores chiral symmetry. 

Fig. 4 illustrates our study of chiral symmetry on the 
original and modified ensembles. As is well-known, xSB 
cannot occur on a finite lattice. Therefore, we measure 
{$ij)){m q ) = (Tr(p + niq)^ 1 ) for a range of quark masses 
m q where finite-size effects are small, and extrapolate 
to m q — > 0. In the original ensemble, (ipip) clearly ex- 
trapolates to a non-zero value which signals x^B. In the 
modified ensemble, the extrapolated value is zero within 
errors: center-vortex removal restores chiral symmetry. 
We expect then the instanton content of the Yang-Mills 
field to be modified also. To check this, we use improved 
cooling jlTj to measure the topological charge of the mod- 
ified field: the striking result is that the removal of center 



vortices always leads to the trivial topological sector. 

We therefore have clear evidence for "center domi- 
nance" : in our modified ensemble, where the center- 
projected field is the trivial vacuum (all links equal to 1), 
the Yang-Mills field shows the IR properties of the trivial 
vacuum, i.e., no confinement, no xSB and no topology. 
The IR properties of the Yang-Mills field appear to be 
determined by its center projection. 

On the other hand, a large number of studies now sup- 
port the alternative scenario of "Abelian dominance." 
We use our approach of center- vortex removal to directly 
assess the relationship between these two scenarios. 

In a first experiment, we construct the Abelian pro- 
jection of our original SU(2) ensemble by gauge-fixing 
to Maximal Abelian Gauge in the usual way then 
identify and remove center vortices from the Abelian sec- 
tor. While the original Abelian-projected ensemble shows 
confinement, with a string tension similar to the non- 
Abelian one, the modified Abelian-projected configura- 
tions do not confine. Therefore, we find no contradiction 
between "Abelian dominance" and "center dominance." 
The latter simply appears more fundamental because of 
the greater reduction of the number of d.o.f. 
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FIG. 5. Size distribution of monopole clusters on the orig- 
inal and modified ensembles. 

In a second experiment, we look at clusters of Abelian 
monopole currents, whose percolation has been identified 
as the signal for confinement Jisj , obtained from the orig- 
inal and modified ensembles. We find that the removal 
of center vortices changes the distribution of monopole 
cluster sizes in a crucial way (see Fig. 5) : whereas in 
the original ensemble, each configuration contains typ- 
ically one very large, percolating monopole cluster and 
many very small ones, the modified ensemble gives a more 
homogeneous size distribution, with a handful of large 
clusters per configuration; these are the remnants of the 
very large one, broken into pieces by the vortex removal. 
Some of them still percolate, even though confinement 
has disappeared. Therefore, we are led to associate con- 
finement with a more specific feature of the monopole 
clusters: monopole current loops which wind around the 
periodic lattice. Such loops can be found frequently on 
the original, confining ensemble, but never on the modi- 
fied, non-confining one. We conclude that: {%) on a finite 
lattice confinement manifests itself in the Abelian sec- 
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tor by the presence of monopole current loops with non- 
trivial topology; (ii) center-vortex removal, which de- 
stroys confinement, always finds the "weak links" of these 
non-trivial loops and breaks them into trivial pieces. 

Now let us consider the issue of gauge-fixing ambigui- 
ties, which was the reason for our initial skepticism about 
the center- vortex idea. These ambiguities come from the 
structure of Q (Eq.[l]) , which has many local maxima, any 
of which can be selected by a local iterative maximization 
algorithm. Each local maximum, or Gribov copy will 
have its own set of P-vortices, differing in number and lo- 
cation. The proposal of p3] is that, no matter which Gri- 
bov copy one chooses, P-vortices are the traces of phys- 
ical center vortices and are roughly located at their cen- 
ter. This argument may account for P-vortices differing 
in location but not in number. To study this question in 
more detail, we magnified the effect of gauge- fixing ambi- 
guities, by stopping our iterative algorithm early, as soon 
as e (Eq.[|) < 10 2 . Thus we not only explore a different 
basin of attraction of Q, but we do not even stop at a local 
maximum. One effect of this partial gauge-fixing is ex- 
pected: the density of P-vortices increases from p « 5.5% 
to w 7.4%, i.e. shows an increase 8p f» 1.9%. The string 
tension measured in the ^-projected ensemble increases 
accordingly: whereas the Zi string tension after "com- 
plete" gauge-fixing (aa 2 ~ 0.075) is a little larger than 
but compatible with the non-Abelian string tension 
(0.0708(11) @, see Fig.2), it jumps to w 0.12 after par- 
tial gauge-fixing. What is remarkable is that this increase 
8(J ~ 0.045 is similar to that obtained by placing the sur- 
plus bp of P-vortices at random, uncorrelated locations: 
8a w — ln(l — 28 p). This makes it plausible that the 
center projection always captures the core d.o.f. relevant 
for IR properties, plus a varying amount of unrelated 
noise ppf] . Indeed, it has been argued th at g auge-fixing 
is not even necessary for center projection [ |19| . In our de- 
scription, center-gauge-fixing acts as a UV noise-filtering 
device, with different Gribov copies letting through dif- 
ferent noise components. Further evidence for this is ob- 
tained by removing from the original SU (2) ensemble the 
center vortices identified after partial gauge-fixing only. 
Just as for "complete" gauge-fixing, we observe that con- 
finement is lost, chiral symmetry restored, and the topol- 
ogy trivial. The only difference is that the modified en- 
semble now has much more short-range disorder. 

In conclusion, we have shown that removal of center- 
vortices from an SU(2) Yang-Mills ensemble causes the 
loss of confinement and the restoration of chiral symme- 
try. One may ask about the connection of the modified 
ensemble {£/'} to the physics of the original SU(2) the- 
ory. Note that only plaquettes of {£/'} at the locations 
of P-vortices differ from those of {U}: hence, as a — > 0, 
their proportion goes to zero as a 2 , since the density of 
P-vortices is physical |T^,[l4|]. Therefore, rewriting Eq.|| 
as U^{x) = Zfj,(x) x Ujjx), we see that the original field 
{[/} has been factorized into a (maximally) central part 



{Z} and a quotient, {£/'}, whose field strength differs 
from the original one only on defects of codimension 2. 
Nevertheless, this small difference alters the physics dra- 
matically: {£/'} has perturbative properties, so that all 
the non-perturbative, IR physics must be carried by {Z}, 
which by definition encodes the center vortices. It would 
be desirable, of course, to formulate an effective action for 
the center-projected theory. Ref. plj] considers an exten- 
sion of the Nambu-Goto action, where the fundamental 
d.o.f. are the 2-dimcnsional random surfaces dual to the 
P-vortices. Ref. Q instead proposes to consider center 
monopoles and their world-lines. We suggest identifying 
a "minimum spanning tree" of negative Zi links respon- 
sible for the P-vortices: perhaps only a subset of them 
form the essential d.o.f. governing the IR properties. 

Finally, our vortex-removal procedure can be used to 
study properties of non-confining non-Abelian fields and 
effects of center-symmetry breaking. For instance, re- 
moving time-like center disorder only would be similar 
to raising the temperature above T c . 

We thank A. Di Giacomo, M. Golterman, J. Greensite, 
T. Kovacs, C. Lang and O. Miyamura for discussions. 
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This is what seems to happen for a generic starting point 
on the gauge orbit. There also exist pathological starting 
points, notably Landau gauge, from which one obtains 
very few center vortices, not associated with confinement. 
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